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Abstract 

We consider the phenomenon of Bose-Einstein condensation in a random growing directed net- 
work. The network grows by the addition of vertices and edges. At each time step the network 
gains a vertex with probabilty p and an edge with probability 1 — p. The new vertex has a fitness 
(a, b) with probability f(a,b). A vertex with fitness (a, b), in-degree i and out-degree j gains a 
new incoming edge with rate a(i + 1) and an outgoing edge with rate b{j + 1). The Bose-Einstein 
condensation occurs as a function of fitness distribution f(a, b). 

PACS numbers: 02.50.cw, 05.40.-a, 89.75Hc. 
Keywords: Random networks 
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I. INTRODUCTION 

Recently there has been much interest in random growing networks |jl 0, [J, both from 
the point of view of theoretical modelling, as well as the empirical study of real networks. 
There is considerable evidence that while traditional Erdos-Renyi random graphs have a 
Poisson degree distribution 4^, real graphs and random growing graphs have a power-law 
degree distribution. 

Of particular interest are directed networks, which can be used to model systems in which 
directed flow takes place. Such networks include the world wide web |5l , the phone-call graph 

n o 

p and the citation graph |7(. 

As a complement to computer simulations and exact solutions of simplified systems, 
thermodynamic formulations have been used to study a number of complex co-operative 
systems including granular media, econophysics, breaking phenomena and many others. In 
a recent paper |8], Bianconi and Barabasi (BB) mapped the different behaviour of random 
growing networks with fitness onto the thermodynamically distinct phases of a free Bose 
gas. The fitness model predicts that, in the large network limit, the fittest node will have 
the most links. This is called the fit-get-rich (FGR) phase. Unlike the scale-free model 
(SFM), in which the degree distribution of the network is power-law, the FGR behaviour 
(phase in the thermodynamic language) has nodes with a very large degree which dominate 
the network. Another phase revealed in j| is the Bose-Einstein condensation (BEC), when 
the fittest node acquires a finite fraction of the total degree. In contrast to BEC, in the 
FGR phase the richest node is not an absolute winner, since its share of the links (i.e. the 
ratio of its degree and the total degree of the network) decays to zero for large system sizes, 
whereas in BEC the winner maintains its share irrespective of system size. 

In this case, the fact that this ratio is constant corresponds to the extensivity property of 
a Bose gas, when the gas keeps a finite fraction of its particles in the ground state. In this 
paper we examine the phenomenon of Bose-Einstein condensation, which was considered for 
an undirected graph in 0,0] on the random directed growing network introduced in jlfl| . 

To do this we work with the model introduced in lpj . introducing a dependence between 
the vertex fitness and an energy. As each vertex has two fitnesses (for in-degree and out- 
degree) then it is necessary to assign two different energy levels to represent each vertex. 
This is detailed in the first sub-section, together with the description of edges in terms of 
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particles assigned to two energy levels. 

The necessary coexistence of two sub-systems forming the network and the formulation 
of an equilibrium condition in a formal way is discussed in detail in sections 3 and 4. 

The occurrence of Bose-Einstein condensation, and its interpretation in a physical frame- 
work makes it possible to describe the directed network in terms of canonical concepts of 
statistical physics, such as a phase diagram and first and second order phase transitions. 
These are introduced in the following sections. 

In the conclusions we emphasize the usefulness of this formulation in providing a clear and 
concise interpretation of different phenomena in networks. So, the possibility of coexistence 
of phases in which the in-degree distribution exhibits a clear winner whereas the out-degree 
shows scale free behaviour, or the simultaneous existence of two different winners (bipolarity) 
and the conditions for its existence are predicted in thermodynamic language. 



II. THE MODEL 

n 

In BB [8], a correspondence between fitness rj and energy e given by 

t] = e-P £ (1) 
was introduced. In our model the vertices have, in general, different fitnesses a and b for 
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the addition of an in- or out-degree. As in 1101, here the model consists of the addition of 
bare vertices (i.e. without edges, but with fitness a for in-degree and b for out-degree) to 
the network with probability p, and the creation of directed edges between vertices with 
probability q = 1 — p. 

A kinetic equation describing the process of directed networks must include the kinetics 
of both in-degree and out-degree simultaneously since they coexist and influence each other. 



Then, the 
j is, as in 



dnetic equation for the mean number of vertices with in-degree i and out-degree 

3, 



dNt ^ h) = jftt^-ij - (i + W;l + W^' 1 ~ {j + 1)Ay + **°W(«> h) - (2) 
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where the fitnesses (a, b) are chosen from the fitness distribution f(a, b). The normalisation 
factors Mi and M 2 are given by 

M i = E(* + i) aAr uM) (3) 

ijab 

and 

M 2 = J2U + l)bN t ,(a,b). (4) 

ijab 

The first term in the first square brackets of Eq.(2) represents the increase in Nij when 
vertices with in-degree i — 1 and out-degree j gain an incoming edge and the second term 
represents the corresponding loss. The second square brackets contain the analogous terms 
for outgoing edges and the last term ensures the continuous addition of new vertices with 
fitnesses a, b. 

The translation of this problem into an energetic formulation is straightforward. We 
associate the addition of a vertex with fitnesses (a, b) to the network with the creation of two 
energy levels; one representing the fitness of the incoming edge and the other for the outgoing 
edge. This means that the edge can be mapped into two separate isolated sub-systems. Then 
the creation of a directed edge corresponds to the creation of two particles, one in each sub- 
system, simultaneously. The particles are created in the energy levels corresponding to 
incoming fitness of the receiving vertex and the outgoing fitness of the emitter vertex, as 
Fig. 1 illustrates. 

Then, using the energy variables, rather that the fitness, we have 

a = e- ft£l , (5) 

b = e-^ 2£2 (6) 

and 

9N ^ £2) = ^^[^^-(i+i^^ (7) 
and the normalisations are re-written as 

M 1 = J2^ + ^ 1£l N l ,(e 1 ,e 2 ) (8) 

ijab 



and 



M 2 = ^(j + l)e-^ 2 ^.( £l , £2 ). (9) 

ijab 



As we distinguish between incoming and outgoing edges, and their respective fitnesses, 
when going to the energy representation the edges and their fitnesses must also be distin- 
guished as belonging to different isolated sub-systems of the same system. Let us denote the 
sub-systems as 1 and 2 So, the creation of a vertex in the network implies the simultaneous 
creation of an energy level in each sub-system, i.e. e% and e 2 . 

The creation of an edge joining two vertices implies the simultaneous creation of one 
particle in each sub-system in the energy levels corresponding to the fitnesses of the vertices 
gaining an in- and out-degree, as Fig.l illustrates. 

It should be noted that we are doing a simultaneous description of two isolated sub- 
systems which compose the whole system (the network). These sub-systems do not exchange 
energy or particles, but are correlated in the sense that the creation of an energy level (a 
particle) in one of the sub-systems implies the simultaneous phenomenon in the other. 

By definition of the model, the total number of particles increases linearly with time, so 

E N i3 (e 1 ,e 2 )=pt. (10) 

ij,ei,£2 

Let us define the reduced moments nil and m 2 by 

M 1 (t) = m 1 t, (11) 
M 2 (t) = m 2 t (12) 

Nij(t) = rnjt. (13) 



and introduce Yiij as 



As shall be seen these magnitudes will be useful to calculate the characteristics of the 
network. Then, following the procedure in lfl|. we can show that 



1 g-a&ei 



mi = — E e 

mi 



{i + l)n i:j +p J e /3l£l f(ei,e 2 )deide 2 



and 



m 2 



m 2 



e~ 2 ^(j + l)mj + P J e-^f(ei, e 2 )deide 2 



(14) 



(15) 



which determines nil and m 2 as functions of p, (3i and f3 2 . 
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A. The coexistence of two systems 



In this representation, as we already pointed out, the addition of a vertex means the 
creation of two energy levels, each corresponding to its fitness. When an edge joining two 
vertices is created, we add one particle in each corresponding level. 

The distribution of in- and out-degree then corresponds to the distribution of particles 
in systems 1 and 2 

is the number of vertices with in-degree i, and 

= = r(j + i)r(i + e «"-"') 

is the number of vertices with out-degree j. For large degree their asymptotic expressions 

are 

g . „ r (l+e««i-«)) (lg) 

and 

^.-r (1+e/3(e2 "" 2)) - (19) 

where we have introduced 

pi k = — (20) 

for k = 1 and fc = 2 as chemical potentials. Unlike BB, the chemical potential is introduced 
here as an exact expression, not as an asymptotic one. 

B. "Equilibrium" condition 

The equations for ni\ and m 2 can be transformed by introducing the generating function 

g(x,y,e 1 ,e 2 ) =J2^ l y : 'n ij (e 1 ,e2) (21) 

where 

g(l,l,e 1 ,e 2 )=pf(e 1 ,e 2 ) (22) 
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with 



and 



Then the equations for m fc (fc = 1, 2) become 

p /■ d£ 1 de 2 f(e 1 ,£ 2 ) 



q J (^Pk^k—^k) — \ 

This means that 



1. (25) 



de 1 de 2 f(e 1 ,e 2 ) _ f de 1 de 2 f{e u e 2 ) 



e /9i(ei-/ii) _ ]_ J e /3 2 (£2-M2) _ J 

Unlike BB, this is a new condition generated by the model. Let us denote it as a generalized 
equilibrium condition. A strong equilibrium condition implies that 

si = e 2 (27) 

and 

/ii = £t 2 . (28) 
This condition will be discussed later. In this formalism the occupation probability is, then, 

Me k ,^) = - q eM£k ^ k) _ 1 , (29) 

for k = 1,2, representing the probability for a given edge to belong to a given vertex 
(probability of a particle to belong to an energy level). If, now, the fitness distribution is 
separable then 

f(e 1 ,e 2 ) = f 1 (e 1 )f 2 (e 2 ) (30) 

and 

g / f" kf{ f = 1 (31) 

separately for each system k = 1,2. In each case, is the solution of this equation, which 
depends on the density of states fk{e). 
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C. Bose-Einstein condensation 



To illustrate, let us start with a particular case, when /(e) = Ce e with < e < £ ma x- 
The normalisation C = as in BB. 
Then the integral in Eq.(25), 

T/ n ^ p(0 + l) f dee 6 

But this integral reaches its maximum value when \x = 0. Then if for a given value of (5 
we have 

M < 1. P3) 
then we have a Bose-Einstein Condensation (BEC). If we calculate the temperature for BEC 
in this particular case, we obtain 

i 

T BE = — \^r(e + 2)((e + i)\ 6+1 ■ (34) 

Unlike BB,the factor 2 appears here. If p = q we recover BB result. Even in this simple 
example, Tbe m &y be different for the sub-systems 1 and 2, if the value of 9 is different. In 
a more general form, we can state the condition for Bose Einstein condensation. Then, as 

t(r „ \ p [ ds 1 de 2 f{e 1 ,e 2 ) 

I((3 k , l*k) = -J eM£k ^ k) _ I = 1 (35) 
for k — 1. 2., if I(Pk, Hk) < 1 then BEC occurs. The conditions for BEC depend on f(ei, £2)- 



D. Phase diagram 

Below Tbe there is BEC. Above it the phase is fit-get-richer as in but here both 
phases must be considered for both systems. Then, a phase diagram can be drawn with the 
temperatures T\ =4^ and T 2 = characterizing each sub-system (Fig. 2). 

Imposing a given dependence of fitness with energy so that 

r de 1 de 2 f(e 1 ,e 2 ) f de 1 de 2 f(s 1 ,e 2 ) 

means that the equilibrium is broken, (i.e. in a computer simulation we would impose 
temperatures that do not satisfy Eq. I2H1)- Then there can be a winner in one of the phases. 
One of the vertices accumulates, say, most of the in-degrees, coming from the other vertices. 
If both phases are in BEC there is a vertex with most of the incoming and another with 
most of the outgoing edges. Then both vertices will be highly correlated. 
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E. Out of equilibrium first order phase transition 

Though a more complicated situations are possible, let us illustrate the case T\be = T 2 be- 
The point A represents that both sub-systems are in FGR phase. If, for example, we make 
the simulations with the parameters given by Al then will be an FGR and a BEC phase 
coexisting, but not in equilibrium since one of them satisfies Eq. EH1 but not the other. 
There is a coexistence of phases out of equilibrium. We might call this an out of equilibrium 
first order phase transition. The motion along AA1A2... leads to the region in which BEC 
occurs for both phases. Here is hard to speak about equilibrium unless A belongs to the line 
T\ = T 2 . All we know in BEC region is that none of the sub-systems satisfies Eq. IHB1 

The concept of equilibrium in this model is purely formal since the sub-systems are not in 
contact, but it is a comfortable tool to make a map of directed networks into Bose-Einstein 
statistics. 



The scale free model 



If one of the systems is such that /(e) = S(e), then Eq. (}3*T]) gives 

T = 1 (37) 
q e-fo - 1 

then 

g t ~ r< 1+e ^ ~ rP+V (38) 

with a corresponding expression for hj. If p = q the scale free model is recovered with a 
degree distribution Pk ~ k~ 3 . Thus the model with directed edges gives the possibility of a 
scale free phase for both in the in-degree and the out-degree. 



G. Second order phase transition 

If in Fig. 3 we move in such a way as to cross into the BEC region through the point 
B then BEC occurs simultaneously for both phases. No coexistence of phases will occur, 
the system goes from FGR to BEC or vice-versa instantly. This is a second order phase 
transition since the system as a whole changes its phase. Then the point B is the critical 
point for the transition. 
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BEC is itself a second order phase transition since it effects the system as a whole, so 
that when sub-system 1 or 2 experiences BEC, the sub-system itself suffers a second order 
phase transition, but our system contains both sub-systems 1 and 2 so that when they are 
in different phases we consider our whole system as composed of two phases. 



III. DISCUSSION 



We have studied Bose-Einstein condensation in a random growing network model. In this 
model, with a formulation considering incoming and outgoing edges, a formal analogy can 
be adopted to describe the change in the behaviour of the directed edges in terms of Bose 
statistics. 

In this formulation, as in that of BB, the temperature plays the role of a dummy variable. 
Real systems are characterised by the functional form of the distribution function of fitness 
and there is no j3 to consider. It only emerges when we translate from the language of 
fitness to that of energies. But the formulation of networks in terms of quantum statistics 
and the introduction of temperature is an elegant and simple description of the behaviour 
of networks under different conditions. In this sense, the temperature plays its role in the 
simulations of the networks, determining the strength of the dependence between the fitness 
and its associated energy, and this dependence is reflected in the behaviour of the network 
by determining its position in the phase diagram. 

Though phase diagrams and phase transitions are not exclusive to thermodynamic for- 
mulations (see, e.g. Dorogovtsev and Mendes ll|), such a formulation is a very powerful 
tool to describe different phases of random directed networks, where the presence of directed 
edges make the description in terms of fitness difficult. Considering this network as a ther- 
modynamic system, makes it simple to understand, provided we appropriately interpret the 
thermodynamic parameters. 
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V. FIGURE CAPTIONS 



1. Correspondence between fitness and energy for directed networks. 

2. Phase diagram. 

3. First and second order phase transitions. 
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